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Abstract 


Explorer 34 plasma data are presented which provide tests of 
current theories of the nature of the boundary between the earth* s 
magnetosheath and the interplanetary medium. It is shown that a 
discontinuous transition between the magnetosheath and the interplanetary 
medium is seen simultaneously in the proton spectrum, the flow direction # 
and the magnetic field intens ity, during both geomagnetically quiet and 
disturbed times. A few boundary crossings were observed for which the 
plasma spectra appear to indicate a diffuse transition; but a detailed 
study of an extreme example of such a crossing, based on simultaneous 
measurements of the spectra, the flow direction, and the magnetic field 
intensity, shows that the apparent diffuseness is actually the result 
of a complex microstructure which consists of a number of closely 
spaced discontinuities that could be due to multiple crossings of a 
bow shock. Observations of the velocity field in the neighborhood of 
the outer boundary of the magnetosheath are presented and are compared 
with results of the fluid theory, which predicts a standing shock. It is 
shown that a along the boundary on the flank of the magnetosheath, the 
ratios of the magnetosheath flow speed to the interplanetary flow speed 
are in close quantitative agreement with a prediction of the fluid 

theory. This prediction is essentially independent of the Mach number and the 
adiabatic exponent. We observe a highly directional flow in the region 
30° < 0 < 100°, where 0 is the sun-earth-satellite angle. The direction 
of this flow is in good agreement with the bow shock theory. 



I* Introduction 


Axford (1962) predicted that a standing shock wave should exist 
permanently on the sunward side of the earth as a result of the flow 
of the supersonic interplanetary plasma past the earth* s magnetosphere. 
Kellogg (1962) independently proposed the same idea and calculated 
the shape of the shock boundary using classical aerodynamic theory 
for supersonic flow past an obstacle* More refined calculations 
of the shock shape were later published by Spreiter and Jones (1963) 
who used the fluid theory with more realistic boundary conditions* 

An alternative theory for the interaction between the interplanetary 
plasma and the earth* s magnetic field was subsequently published 
by Bernstein et al (1964), Fredericks et al (1965), and by Scarf et 
al (1965)* This free-flow theory implies a gradual, broad, disordered 
transition between the magnetosphere and the interplanetary medium. 

Ness (1967) in his review of experimental work on the solar wind 
geomagnetic field interaction has noted that there is, at present, 
disagreement as to which of these two theories is applicable. 

There are a number of observations which support the shock theory; 
Boundary thickness . The magnetic field data of Ness et al (1964) and 
Heppner et al (1967) consistently showed a thin boundary between the 
magnetosheath and the interplanetary medium* Heppner et al (1967) 
have determined that the thickness of the boundary seen in the magnetic 
field is less than ^250 km* Argo et al (1966) reported that they saw 
no evidence for a diffuse boundary in the plasma data from Vela 2. 
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Boundary shape 0 Dryer and Heckman (1967) have reviewed both plasma 
and magnetic field observations of the position of the boundary and 
have found that the observations are consistent with the shock theory. 
Velocity field o Examining the flow near the boundary at about 60° and 
300° solar ecliptic longitude and between + 56° ecliptic latitude, 

Argo et al (1967) found that for 13 crossings the ratio of the magneto- 
sheath flow speed to the interplanetary flow speed agreed with the 
theory pf Spreiter et al (1966) and Dryer et al (1966). They observed 
an ordered, directional flow, as the shock theory predicts . 

Temperature ratios ** Spreiter et al (1968) have noted that the observations 
of the temperature ratios reported by Argo et al (1967) agree with 
the shock theory* 

Observational evidence against the shock theory and in support 
of the free-flow theory was presented by Wolfe et al (1966a, 1966b): 
Boundary thickness <, Wolfe et al ( 1966b) reported that they observed 
a diffuse boundary for more than half of the first twenty-one orbits 
of IMP 1. They suggested that the solar wind may be only marginally 
capable of producing a detached bow shock and that the shock may not 
form during geomagnetically quiet periods when the solar wind flux 
and speed is low* 

Velocity field , Using data from four satellites (two near the subsolar 
point and one far out on each flank of the magnetosheath, near the 
magnetopause), Wolfe et al (1966a) found no increase in the bulk speed 
as the plasma flows around the magnetosphere. This is inconsistent 
with the shock theory. 
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Clearly, the evidence in support of the fluid theory is substantial 
and convincingo The evidence against the fluid theory concerns the 
boundary thickness and the flow velocity behind the boundary. In 
this paper we examine the problem of the boundary thickness, using 
simultaneous plasma and magnetic field data, and we examine the plasma 
velocity taking care to consider several of the variables involved. 

Our results provide additional support for the fluid theory. 



4 


II . Instrument, Orbit and Data Reduction 

Explorer 34, launched on 27 May 1967, carries a plasma 
experiment provided by the present authors. Since a full technical 
description will be presented elsewhere, (Ogilvie, Mcllwraith and 
Wilkerson, 1968) only a condensed version will be given here. 

Instrument * The experiment makes observations of the spectrum of ions 
in the energy per unit charge range 310 eV per charge to 3100 eV per 
charge. Fifteen energy intervals of width 4.57 f a are sampled, and their 
center energies form a geometric series with a ratio of 1.24. An 
observation consists of counting ions whose energies fall in the pass 
band, and whose trajectories within the solid angle, of the analyser 
during one satellite revolution. As the sensitive angular cone of the 
experiment crosses the anti-solar direction, the energy sampled is 
changed to the next in order. The experimental program is shown in 
Figure 1, where it will be seen that in the course of nine telemetry 
sequences or about three minutes, three spectra are taken. The first 
of these is the energy spectrum of protons. The second is the energy 
spectrum of doubly charged helium ions and the third consists of 15 
successive readings on the same proton energy 'step’; helium data 
will not be discussed in this paper. 

The time resolution of the experiment is specified 
as follows: A single observation of the ions of one species at one 

energy per unit charge setting occupies 2.6 seconds, a complete 
spectrum of ions of a single species requires approximately one minute 
to record and these are repeated every 3 minutes. This program continues 
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without interruption. The ion to be studied, H + or H e ++ , is selected 
unambiguously by setting the velocity selector interposed between the 
analyser and the detector. 

As an example of the operation of this device, Figure 2 shows 
a typical solar wind proton spectrum. The degree of separation achieved 
between the two ions may be expressed as the ratio of the number of 
counts with the instrument adjusted to observe ion A, per incident ion 
B, values of about 1CT^ are typical. The background seen in Figure 2 
is a function of the operation of the computer which compresses the 
data before transmission, and is not caused by the detector. It is 
completely explained and determined during the recording of each 
spectrum so that it may readily be substracted before data reduction; 
there has been no change in this rate during the first 78 days in orbit. 
The threshold indicated in this figure represents a level of three 
standard deviations above the background level. Counting rates exceeding 
the threshold are regarded as significant. 

The sensitive cone of the instrument, which is two degrees 
in azimuth and +9° in a plane containing the spin axis, rotates with 
its axis of symmetry in the plane of the ecliptic. By recording the 
number of particles detected by the instrument in successive lb^* 13 of 
a revolution, a series of 16 numbers, C^, is obtained. Counting 
starts when the optical aspect sensor, looking out into a cone of 
angle 1° whose axis is normal to the face of the sector opposite to 
that occupied by the plasma detector, sees the sun. The 16 subtotal 
counts then describe the azimuthal distribution of flux recorded as 
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a histogram of 16 bars. One of the telemetered quantities is the 
number of the angular sector from which the greatest subtotal count 
C ? ax was recorded. As seen from Figure 5, and taking into account 
the 15^ degree displacement of the sensitive cone of the instrument 
and the 4 degree aberration correction, cold plasma flowing radially 
from the sun will be "seen” by the instrument in sector 9, 19^ degrees 
from the boundary between sectors 8 and 9. Thus plasma flowing from 
a direction displaced by more than about 3 degrees westward of the 
sun will be detected in sector 10, later in the rotation. A deflection 
of the flow direction towards the east by more than 19^ degrees will 
cause detection during the traversal of sector 8. These qualitative 
angular estimates can be refined if the shape of the angular distri- 
bution is assumed; this point will not be pursued here. 

Orbit. Figure 4a shows the trajectory of Explorer 34 during orbit 2, 
and it shows the position of apogee on orbits 2 to 20. It can be seen 
that the satellite was launched at a sun-earth-satellite angle 0 S g=lO6 u 
and apogee moved at the rate of approximately 4.3° per orbit toward 
the subsolar point , 0 SE decreasing. The orbital period is approximately 
4,3 days. Apogee and perigee are 34 Rg and < 1 R £ , respectively. 

The orbital plane is nearly perpendicular to the ecliptic plane. 
Projections of the trajectory in the X, Z plane are shown in Figure 4b 
for several orbits; this is the projection which is seen when viewing 
away from the sun along the earth-sun line. 

The spin axis of the satellite is oriented 2.5° to the normal of 
the ecliptic plane; in what follows, we shall neglect this small angle. 
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Reduction of Data to Obtain Fluid Parameters Details of the data 

reduction program are given in a report by Ogilvie et al.(1967). Briefly, 
the fluid parameters n, u and T are obtained from the observed spectra as 
follows (see Figure 2). 1) The working spectrum consisting of NOC adjacent 

bars is obtained by setting counts equal to zero in those energy 
channels in which the counts are less than some threshold. 2) The 
remaining counts are corrected for background and the counts in 
channel i (i = LF, NOC-1) are then converted to a value ot the distri- 
bution function (dn/dv) by means of the equation 

2 

(— ) = constant x (C-j/V,-) 
qv x 1 

which is obtained from the solution of the unfolding problem that is 
given by Ogilvie et al. 3) The first three bars in the new working 
spectrum are fitted by the function 

f(v) = vexpLa^+3-^v+v^v ] (1) 

to obtain a representation to the actual observed distribution 
f (v) = ^ (v) between channels LF and LF+1; the next three bars 
(LF+1 , LF+2, and LF+3) are similarly fitted to obtain an approximation 
to f(v) between LF+1 and LF+2, and this process continues until an 
analytical approximation to f(v) is obtained by fitting all of the 
observed channels to segments of the function given by (1). 4) The 
fluid parameters are calculated by integrating the function f a (v) 
which is obtained by the fitting procedure: 
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n = J q f a (v) dv (2) 

<V> = i v f a (v) dv (3) 

T = fn { i v2 fa(V ^ dv " n <v>2 } (4) 

This procedure gives reliable values of the fluid parameters in the 
interplanetary medium when the number of bars in the working spectrum 
is greater than or equal to three. In the magne to sheath , however, 
the parameters obtained in this way must be used with caution. The 
reason is that in the magnetosheath, where the temperature may be 
high and the bulk speed low, the angular spread of the distribution 
function may be greater than the aperature of the instrument. Con- 
sequently, the detector may not accept all of the particles, and the 
observed distribution f a (v) (which is an integral of f (v) over a 
fraction of velocity space) may be significantly different from the 
actual speed distribution function f(v) (which is an integral of f(v) 
over all of the velocity space). This effect has been studied in 
detail for a convected maxwellian distribution. Equations (A 18) and 
(A 26) in the Appendix give the values of dn/dv which should be observed 
for a given density, temperature, and bulk speed if the flow direction 
makes an angle 9 Q with the axis of the acceptance cone of the detector 
which has a half-width A0 ^ 9°. The n c , T 0 , and u Q which should be 
observed for a given n, % and u have been computed from (3), (4), and 
(5) by setting f a = dn/dv where dn/dv is given by (A 18) and (A 20) . 
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The results are shown in Figure 5. ]t can be seen that the mean 
speed <v> is essentially independent of 0 Q . The mean speed is generally 
higher than the flow speed v, but this difference is less than 5% when 
the temperature is less than 5 x 10^°K. The density 
is relatively insensitive to 0 O when T ^ 5 x 1C)5 0 k, but it is 
more than a factor of two lower than the actual density when T ^ 10^°K. 
Thus, in the magnetosheath we can set <v> = u when T ^ 10^°K with an 
error less than 10%, but the errors iu density may be large if no 
corrections are made. In the interplanetary medium the fluid parameters 
are correct to < 10% if T < 5 x 10^°K. 
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III . Existence of a Discontinuous Transition Between the Magnetosheath 
and the Interplanetary Medium . 

The bow shock theories which were mentioned in the Intro- 
duction imply that the boundary between the magnetosheath and the 
interplanetary medium should be very thin. Hydromagnet ic shock theories 
such as that of Kennel and Sagdeev (1967) give a shock thickness some- 
what larger than an ion cyclotron radius, viz. L ^ few x 100 km ^ .1 . 

Shock theories involving ion wave instabilities, 

such as that of Tidman (1967), give a shock thickness on the order of 
the bulk-speed divided by the ion plasma frequency in the interplanetary 
medium, viz. L % few x 10 km ^ .01 . The magnetic field measurements 

of Heppner et al. (1967) show that the major change during a bow shock 
crossing occurs within 30 seconds, from which they conclude that the 
thickness of the bow shock is less than 1 . Thus, although there 

is no agreement as to precisely what the shock thickness should be and 
although it is difficult to measure the thickness because of motions 
of the shock, proponents of any shock theory would probably agree that 
an upper limit to the apparent shock thickness is % ,1 R £ . Since the 
Explorer 34 satellite typically moves approximately 1 R^ in several 
minutes near the outer boundary, and since the plasma analyzer records 
successive proton spectra at three minute intervals, the shock theory 
implies that a shock crossing viewed on a scale of several hours would 
appear as a discontinuous transition between the interplanetary medium 
and the magnetosheath in the overwhelming majority of cases. On the 
other hand, the free flow theory apparently implies a very diffuse 
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outer boundary with a scale on the order of a few earth radii (see 
Wolfe et a 1 , 1966b) . 

The purpose of this section is to show that we observe a 
discontinuous transition between the magnetosheath and interplanetary 
medium during both quiet and disturbed times. We shall show that 
this transition is seen simultaneously (within the three minute reso- 
lution of the plasma probe) by the magnetometer and the plasma probe. 

Quiet Time Transition . We shall first describe a transition between 
the magnetosheath and the interplanetary medium which was observed when 
the satellite moved from 17.4 Rg to 14.0 at an earth-sun-apogee angle 
of ~58° in the interval 1600 UT-1900 UT on July 10, 1967. We wish to 

emphasize that July 10 was a very quiet day. Its IAGA designation is 

QQ , the three-hour Kp indices for this day are 0 + , 0 + , 1, l", 0 + , 1, 2. 

The lowest graph in Figure 6 shows the magnetic field in- 
tensity as a function of time. The low values of B in the interval 

1600 UT - 1804 UT are characteristic of the quiet time interplanetary 

magnetic fields, and the large magnetic field intensities in the 
interval 1804 UT - 1900 UT are characteristic of the outer magneto- 
sheath (see Ness, 1967). The detailed magnetic field data show that 
. the transition between these two regions occurred in less than 20 
seconds. During this time the satellite moved .005 R £ . 

The plasma flow direction in the interval 1600 UT - 1900 UT 
on July 10 is shown above the magnetic field data in Figure 6. Prior 
to 1804 UT the flow is into channel 9 as one would expect from the 
nearly radial solar wind. After J 804 UT the flow is into channel 11. 
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Section V shows that such a flow direction is consistent 

with that predicted by the theory of Spreiter et al. (1966) 

for the region of the magnetosheath just behind th£ bow shock at 0gE = ^8° • 

Thus, the change in flow directions in Figure 6 indicates a discontinuity 

at 1804 UT. The width of the transition determined in this way is 

less than three minutes which corresponds to L < .05 R E . 

Plasma spectral data for 1600-1900 UT on July 10 are summarized 
in the top graph of Figure 6, which we refer to as a plasma spectrogram 
or "plasmagram". The vertical axis of the plasmagram shows speeds 
corresponding to the centers of the velocity channels of the plasma 
detector. Each column of numbers in the plasmagram shows a measure 
of the counts (see Table 1) which were recorded in the various energy 
channels during a one minute interval. In other words, each column 
represents a complete proton spectrum. The plasmagram is a series of 
these spectra ordered on a linear time scale. The narrow spectra in 
the interval 1600 UT - 1804 UT are characteristic of the interplanetary 
spectra which are observed during quiet times. The low bulk speed 
(u ps 300 km/sec) agrees with the quiet time interplanetary bulk speed 
reported by Wolfe et a 1. (1966b) and by Hundhausen et al- (1967) and 
also with theoretical predictions. The low temperature (T J> 2 x 10 4 °K) 
is also typical of the quiet interplanetary medium. The broad spectra 
after 1804 UT are undoubtedly magnetosheath spectra. We do not observe 
the low energy ends of these spectra because the peaks occurred near 
the lowest energy step of the instrument, so we cannot compute a mean 
speed or temperature. Nevertheless it is clear that the mean speed 
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has decreased and the temperature increased considerably after 1804 UT. 
These spectra show a variability which is characteristic of magneto- 
sheath spectra which we observed on other orbits. 

In summary, the combined plasma and magnetic field data in 
Figure 6 show a d iscont inuous transition (boundary thickness < . 1 Re) 
between the interplanetary medium and the earth’s magnetosheath during 
a geomagnetic quiet time. This agrees with the fluid theory. It is 
clearly incompatible with the free flow theory of the tran- 
sition region. It is also inconsistent with the hypothesis that 
a discontinuous transition exists only at disturbed times. 
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A "Diffuse" Transition During the first 15 orbits of Explorer 
34 many discontinuous transitions were observed as the satellite 
moved between the interplanetary medium and the magnetosheattu 
However, there are a few cases (~ 5 %) in which the plasma spectral 
data alone show a transition between the magnetosheath and the 
interplanetary medium which appears to be diffuse „ An extreme 
and unusual example of such a transition is shown in the plasmagram 
in Figure 7. 

In Figure 7 the broad, fluctuating spectra with low counts, 
which were observed during the interval 2110-220 UT, are characteristic 
magnetosheath spectra 6 The narrow spectra with higher counting 
rates, which were observed between 2230 and 2330 UT, are characteristic 
interplanetary spectra. There appears to be a gradual transition 
between these two types of spectra in the interval 2200-2230 UT„ 

During this time the satellite moved 0 J This apparently diffuse 

boundary occurred during a geomagnetically quiet period in which 
the 3 hour Kp index remained near 1 for 18 hours. Thus, in the 
absence of magnetic field data and without information about the 
plasma flow direction, it is easy to see how one might infer, with 
some justification, that a diffuse outer boundary might exist during 
a geomagnetic quiet periodo However, we shall now show that more 
comprehensive data for this period indicate such a conclusion is 
probably incorrect., 
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The magnetic field intensity shown in Figure 7 indicates 
that the satellite was in the magnetosheath between 2100 UT and 
2217 UTo At 2217 UT the magnetic field data show a discontinuous 
decrease to typical interplanetary intensities. This discontinuity 
was accompanied by a discontinuity in the flow direction from 
channel 11, which is typical of the magnetosheath, to channel 10, 
which is commonly observed in the interplanetary medium. A change 
in the character of the spectrum can also be observed near 1217 UT. 
The spectrum change alone is ambiguous, but the combination of 
observations just described definitely implies the crossing of a 
thin boundary. Another discontinuity, which is the complement 
of that just described, can be seen at 2323 UT. .Here the magnetic 
field increases, the direction changes from sector 9 to 10, and the 
spectrum broadens, indicating an abrupt transition from the 
interplanetary medium to the magnetosheath. 

In the interval between 2231 UT and 2234 UT (marked M in 
Figure 7) the magnetic field data are characteristic of neither the 
magnetosheath nor the interplanetary medium. This observation is 
not inconsistent with the shock theory, for it is possible that the 
shock moved rapidly in and out during this interval so that both 
magnetosheath and interplanetary spectra were recorded but not 
clearly time resolved. Only one spectrum was recorded during this 
interval and the counts were too low to draw conclusions from it. 
Another "mixed" region can be seen between 2330 UT and 2335 UT in 
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Figure 7. Here two spectra can be seen, but again the counts 
are too low to determine whether they indicate interplanetary 
or magnetosheath conditions* However, the magnetic field data 
strongly suggest that both conditions were observed in this 
interval as a result of multiple crossings. 

Finally, let us consider the peculiar spectrum at 2358 UT„ 

The spectrum is very broad which is indicative of the magnetosheath, 
but the flow is from channel 9, which is indicative of the 
interplanetary medium. This peculiar dilemma is readily resolved 
by referring to the magnetic field data. The magnetic field 
intensity shows that a discontinuous transition from the magneto- 
sheath to the interplanetary medium occurred while the spectrum 
was being recorded. This implies that the spectrum should be a 
mixture of a magnetosheath spectrum and an interplanetary spectrum, 
consisting of low counts which were recorded in the low energy channels 
before the boundary crossing and higher counts which were recorded in 
the higher energy channels after the crossing. The observed spectrum does 
show these characteristics. Moreover , • the observation of flow 
from channel 9 is consistent with t \ts interpretation, for we 
plotted the flow direction which was recorded when the highest 
count was measured, and in this case the highest count was measured 
in the interplanetary part of the mixed spectrum. 

We conclude that although the unusual plasmagram shown in 
Figure 7 may suggest, at first sight, a diffuse transition from 
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the magnetosheath to the interplanetary medium, the apparent 
diffuseness actually is the result of a complex microstructure, 
which is revealed by simultaneously examining the plasmagrai?, 
the flow direction, and the magnetic field intensity. The 
microstructure consists of a series of discontinuities, which we 
interpret as multiple crossings of a thin boundary between the 
magnetosheath and the interplanetary medium. 
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Typical Transitions On May 30, 1967, an unusual combination* of 
"typical" bow shock crossings was observed in a six hour interval. 

This is shown in Figure 8, which we shall now discuss in detail. 

At 1429 UT on May 30, 1967, a sudden commencement was observed 
at magnetic observatories around the world. At this time (to the re- 
solution of our plasma analyzer), a marked discontinuity in the plasma 
and magnetic field parameters was observed by IMP F which was near 
the bow shock, at R ^ 33.0 R £ , 0g^lOCfl(See Figure 8). Prior to this 
time the geomagnetic field was quiet (for the first four 3-hour 
intervals on May 30, the Kp indices are 2", 2, 2 + and 2"). After 
1429 UT on May 30, the geomagnetic field was distrubed (for the last 
four 3-hour intervals on the date the Kp indices are 6“ , 6, 6 + f 6"). 

Thus, Figure 8 shows data for both a quiet period and a disturbed 
period. 

During the quiet period in Figure 8, we see two types of re- 
gions. Between 1100 UT and 1308 UT, and between 1336-1406 UT, we 
observe a cool plasma with a low magnetic field moving nearly radially 
from the sun at 415 km/sec. This is clearly a "quiet” interplanetary 
plasma. Between 1308 UT at 1336 UT, and between 1406 UT and 1429 UT, 
the data in Figure 8 indicate that the satellite was in the quiet 
magnetosheath. Note that discontinuous transitions occurred simultaneously 
in the plasma and magnetic field data at 1308 UT, 1336 UT, and 1406 UT. 

We assume that these are the result of motion of the bow shock back 
and forth past the satellite. We see such multiple crossings on 
nearly every orbit. 
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During the disturbed period in Figure 8 (1429 UT - 1700 UT) , 
we again see two types of regions. Between 1429 UT and 1557 UT> and 
between 1609 UT - 1700 UT,we see a hot, high speed plasma with a 
large magnetic field whose velocity fluctuates appreciably. This 
region is clearly the disturbed interplanetary medium. The other 
region, between 1557 UT and 1609 UT in Figure 8, is probably the dis- 
turbed magnetosheath. This is suggested by the very broad spectra, 
the flow direction (sector 11) . the high magnetic field, and the low 
bulk speed. Note that discontinuous transitions between the disturbed 
magnetosheath and the distorted interplanetary medium occurred 1557 UT 
and at 1609 UT. 

The above results lead us to the following interpretations 
of the events in Figure 8. Prior to 1409 UT, the satellite was near 
the bow shock, and the plasma near this region was rather quiescent. 

The bow shock was locally moving in and out, and twice it moved out 
past the satellite so that twice the instruments sampled the magneto- 
sheath. Just before 1429 UT, the satellite was in the magnetosheath. 

At 1429 UT a discontinuity (presumably a shock wave which was propagating 
in the interplanetary medium) arrived at the satellite, and the disturbed 
plasma behind this discontinuity drove the bow shock past the satellite, 
toward the earth. Thus, at 1429 UT the satellite emerged from the quiet 
magnetosheath into a disturbed region in the interplanetary medium. 

The bow shock continued to move in and out. At 1557 UT it moved out 
past the satellite and it returned at 1609 UT, giving us a brief look 
at the disturbed magnetosheath. 
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This event is unusual^but it also serves to illustrate the 
features of the outer boundary which are commonly observed in the 
plasma and magnetic field data, name ly^ mu ltiple , discontinuous tran- 
sitions between the interplanetary medium and the magnetosheath during 
both quiet and disturbed times. 
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IV. Quiet-time Flow Speeds Near the Outer Magnetosphere 

Boundary 

Theory . Spreiter et al (1966) have used their steady state 
aerodynamic model of the earth's bow shock to compute u s /u., 
n g / n i a nd Tg/T^ at positions along the bow shock for quiet 
interplanetary conditions; here u, n and T are the flow speed, 
density and temperature, respectively, and the subscripts s and i 
refer to the magnetosheath and interplanetary medium, respectively. 
The authors show that these ratios depend only on the interplanetary 
free-stream Mach number and the adiabatic exponent y when the 
ratios are considered as functions of X/D, where X is the distance 
measured from the earth toward the sun along the earth-sun line, 
and D is given in units of earth-radii by the equation 


D=^k 


1/3 


2 1/6 

(2tt p t u t ) 


(5) 


where p. is the interplanetary proton density and H is the earth's 
i Po 

mean permanent magnetic field at the geomagnetic equator (Hp Q ~ .312). 

It is found that the density and temperature ratios are rather 

sensitive to the values of Moo and y. Since y is essentially 

unknown and since varies with time, we shall not consider density and 



temperature in this report. On the other hand, it is found that the 
theoretical flow speed ratio is not sensitive to the values 
of M and y in the region X/D <.6. Thus the speed ratio gives 

OO ' 

a means of testing the shock theory which does not require an 
exact knowledge of y 9 and it gives essentially a single curve 
with which we may compare speed ratios that are observed during 
a variety of interplanetary conditions. 

To compare the theoretical speeds with the observed speeds, 
it must be understood that the theory as presented by Spreiter 
et al (1966) gives the flow speed u which is essentially the speed 
with which the distribution function is convected, whereas the 
measurements give the mean speed of the observed particles 
which depends on both the convection and the thermal motion of 
particles « Specifically, the flow velocity u is defined in terms 
of the complete, normalized distribution function by the equation 


u = \ v f(v) d v 


( 6 ) 


where dv indicates that the integration extends over all of velocity 
space o This definition of flow speed is introduced when the field 
equations are derived by taking moments of the Boltzman equation. 

The mean speed, on the other hand, is defined by the equation 

' CO 

< v > = (v f (v)dv, 


( 7 ) 
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where f(v) is the speed distribution which is obtained by 
integrating f(v) over all directions. 

To relate |u| to < v > requires an assumption about f(v) . 
Since the fluid theory assumes a steady, equilibrium flow, the 
proper assumption is that f(v) should be maxwellian. 

Using cartesian coordinates, it is easy to verify that if the 

— *■ 

distribution function f(v) is maxwellian in a frame moving with 
speed u, then the magnitude of the flow velocity calculated from 
(7) is equal to u, i.e., |u| = u. In the Appendix it is shown 
that for a maxwellian distribution which is convected with speed 
u, 

— { - \ 2 

f(v) = ^2 = n<$ e 2kT V U (8; 

dv 2kT u 

(see (A16). Substituting this into (7) we find that < v > is 
related to u by the equation 

< v u [l + % ($ 2 + 0 [ <$ 4 ] ) 

where 

v -yzsr 

m 

It is convenient to write (11) in the form 


(9) 


( 10 ) 


w(km/sec) = 13 

10^ 


(ID 
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Note that the most probable speed v^, which is obtained by 
setting the derivative of f(v) in (8) equal to zero, is equal 
to the mean speed to order (w/u) . Thus the mean speed, the most 
probable speed and the flow speed u are equal when w/u < < 1. For 
example, in the interplanetary medium T ~ 5 x 10^ and u « 400 km/sec, 
so ~ *003 and <v> « ^ |u|. In the magnetosheath, the 

difference between <v> and |u| may be significant. Since the 
higher order terms in (9) depend on the magnetosheath temperature, 
which is highly variable and cannot be accurately determined directly, 
a simple comparison between experiment and theory is not possible 
when the higher order terms are significant. From Figure 15 of 
Spreiter et al (1966) it is apparent that the higher order 
terms will be large when the interplanetary temperature and 
Mach number are high. For the purpose of testing the shock 
theory, we shall consider only bow shock crossings for which < 

5 x 10^°K. In this case, the theory predicts that T g < 10^°K in 
the region X/D < .6 and T g < 5 x 10^°K in the region X/D < -.5. 

2 

For these numbers, we find that the magnitude of the term 

s 

is less than .05 and consequently <v g > ~ u g . 

To summarize, we have noted that the bow shock theory gives 
a single curve for the flow speed ratio in the region X/D < .6, 
which is not sensitive to the values of and assumptions 
concerning y. We have shown that during "quiet times" 

(T < 5 x 10^°K, u < 450 km/sec) the flow speed can be determined 
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with an accuracy of a few percent. Thus, quantitative 
agreement between the measured and predicted flow speed ratios 
in the region X/D < . 6 is a necessary condition for the validity 

of the shock theory. We shall now present plasma data from 
the first ten orbits of Explorer 34 which show that the theory 


does indeed meet this test. 
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Observations The measurements of flow speed rations u s /ui 
for bow shock crossings which we shall now discuss were obtained 
during the first ten orbits of Explorer 34 when the satellite 
moved from 0 « 106° to 0 « 60°. Many boundary crossings were 
observed, since multiple crossings occurred on nearly every orbit, 
but only a few of these crossings are suitable for the test 
described above. We select data according to the following criteria; 

1) The temperature on the interplanetary side of the 

4o 

boundary must be less than 5 x 10 K. 

2) At least five identifiable spectra on each side of the 
boundary must be available. Thus, crossings which occurred 
within 15 minutes of one another are not considered. 

3) The flow speed in the magnetosheath must be sufficiently 
high (u g > 270 km/sec) that the peak of the distribution 
function can clearly be seen. 

When these criteria are applied we are left with boundary 
crossings for the first ten orbits for which we can compute u /u. 

S 1* 

as a function of X/D. For each of these boundary crossings we 
have at least five values of u on both sides of the boundary. 

We have computed the ratio Ug/u^ for each crossing in two ways. 

1). By averaging the five values of u Q and u^. This 
method suffers from the disadvantage that the 
plasma parameters may change during the fifteen minutes 
required to obtain five spectra. 
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2). By using the values of u and u. appropriate to 

S 3- 

two single spectra. In order to avoid problems with 

"mixed 11 spectra we have not used the two spectra 

between which the boundary crossing took place, but 

have employed the spectra preceding and following these. 

We have verified that both methods of computing u s /u^ give 

essentially the same result. This gives a check on the validity 

of the steady flow assumption. 

The results are shown in Figure 9 where Ug/u^ are plotted 

as a function of X/D, where D is computed from n^ and u^ using (5). 

The error bars in Figure 9 are simply estimates which were 

obtained by equating Au to one-half the difference of the maximum 

and the minimum values of the five measurements of u and calculating 

A(u /v.) from the equation 
s l 



% 


( 12 ) 


The largest error bars indicate either a low u g , in which case we 
have set u s equal to the most probable speed, or a fluctuating 
flow speed in the sheath and interplanetary medium. 

The theoretical curve in Figure 9 is taken from Figure 14 of 
Spreiter et aL (1966). It is evident that the observed speed ratios 
are in quantitative agreement with the theory within the errors of 
measurement. The measured and observed speed ratios differ by 
less than 10%. 
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V, Flow Directions Near the Outer Boundary . 

In this section the flow directions which were observed in 
the magnetosheath near the outer boundary in the interval 30° ^ N 100° 
are compared with the corresponding flow directions predicted by the 
shock model of Spreiter et al. (1966) 

In Section II it was explained that the flow direction is 
measured by determining in which of 16,22.5° sectors the maximum number 
of counts was recorded during one satellite rotation. The relation 
between the sector number and the actual direction with respect to the 
sun is shown in Figure 4. 

The flow direction data which will be presented were selected 
according to the following criteria: 

1) The data were obtained within ,5 of a clearly defined 
boundary crossing which was within 2.5 Rg of the theoretical 
boundary for = 8, y = 2. 

2) The flow direction determined for a given rotation is 
used if and only if the total number of counts recorded 
during that rotation is greater than 6550. The choice of 
the number 6550 is based on a study of the onboard computer 
which determines the flow direction. This criterion eliminates 
much of the data in the magnetosheath, where the counting 

rate is low, but it guarantees that the data which are 
selected are accurate. 

Seven of the first 18 orbits of Explorer 34 gave more thar. 
one measurement of the flow direction which satisfies the two criteria 
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above. The measurements are summarized in Table 2. Accurate data 
were obtained on orbits 2, 3, 5, 7, 12, 14, and 18. On each orbit 
the most probable flow direction was into sector 11. However, a few 
observations of flow into sector 10 and sector 12 were recorded, this 
could indicate a small variability of the flow direction or a residual 
instrumental effect. It is useful to define a mean sector for each 
orbit by the equation 

mean sector = E Nj_ S±/ ENj_ (13) 

i 

where is the number of observations of flow into sector S^. For 
example, on orbit 2, 16 measurements of the flow direction were made 
which satisfied the selection criteria; 4 measurements indicated flow 
into sector 10, and 12 measurements indicated flow into sector 11, 
so the mean sector 15 10.7 + .4. Values of the mean sectors, given 
in Table 2, show that the flow was predominantly into channel 11; 
the RMS errors in the mean sectors shows that the flow was highly 
collimated. 

The measured most probable, flow directions are compared with 
the theoretical flow directions of Spreiter et al (1966) in Figure 10. . 
This Figure shows the theoretical boundary and the theoretical stream- 
lines for the model with M» = 8,y = 2,D =10 Rg (Spreiter et al have 
shown that the streamlines for ^ = 8,y = 5/3 and D =10 Rg are 
essentially the same). The measured flow direction for each boundary 
crossing is shown in Figure 10 as flow into a narrow sector which is 
placed at the shock boundary at the appropriate 0. It is evident that 
the predicted and measured flow directions are in good agreement. 



VI. Conclus ions 
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We have presented plasma and magnetic field data from 
Explorer 34 which are relevant to the study of the nature of 
the boundary between the magnetosheath and the interplanetary 
medium. The results may be summarized as follows: 

1) The boundary is characterized by a discontinuous change 
in the proton spectrum, the flow direction and the magnetic 
field intensity. The discontinuities in the plasma and 
magnetic field are simultaneous within the 3-minute 
resolution of the plasma analyzer. The thin boundary is 
observed during both geomagnet ically quiet and disturbed 
times. There is no evidence for a diffuse boundary as 
suggested by the free flow theory. 

2) The flow speed in the region -.65 < . 65 agrees 

quantitatively (within 10%) with a prediction of the 
shock theory of Spreiter et al (1966) which is essentially 
independent of the Mach number and the adiabatic exponent. 

3) The measured flow direction in the region 30° < 0 <; 100° 
agrees closely with the theory of Spreiter et al (1966). The 
flow is highly directional, as predicted. 

Recall that other measurements which have been published - 

% 

plasma measurements of the position of the outer boundary (Dryer 
and Heckman, 1967) , magnetic field measurements of the position of 
the boundary thickness (Heppner et al. 1967), and the plasma measure- 
ments of the change in density, bulk speed, and temperature across 
the boundary at 15 R E (Argo et al. , 1967) - are also in quantitative 
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agreement with predictions of the shock theory. These observations 
together with the results of this paper provide strong, quantitative 
support for the fluid theory of the interaction of the solar wind 
with the earth 1 s magnetosphere. 
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APPENDIX 


The purpose of this appendix is to present the convected 
maxwellian distribution function f(v) in spherical coordinates, 
and to calculate the functional form of the distribution, dn/dv, 
which the IMP F plasma probe records in a convected maxwellian 
plasma. 

In a coordinate system which moves with the plasma the 
maxwellian distribution has the form 


f(v) = 


d 3 n 

dv x dv y dv 2 


= A 


r 2, 2 ^ 2-1 

-OfLv +v +v J 
x y z 


A1 


where 


a = 


2kT 


A2 


and 


A = n ( 


m 


3/2 


2nkT 


A3 


Consider a fixed coordinate system such that the v z axis is 
parallel to the spin axis (approximately perpendicular to the 
ecliptic plane) and the bulk velocity is in the v y” v z plane. In 
a spherical coordinate system the bulk velocity makes an angle 


0 q with the ecliptic plane and a particle has a velocity given by 
(v, 0, 0). The transformation equations which relate ( v x > v y > v z ) 
to (v, 0, 0) are 
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v = vcos 0 sin 0 
x 


v = vcos 0 cos 0 - u cos 0 A4 

y ° 


v„ = vsin 0 - u s in 0 
z o 


It follows that 

v 2 +v, 2 +v 2 = v 2 +u 2 -2uvcos0cos0cos0 - 2uvsin0sin0 r , A5 
x y z o u 

and it can readily be shown that the Jacobian of the transformation 
is 


S(v x ,v v ,v z ) 2 

* = v cos 0 

S(v ,0,0) 


A6 


Substituting A5 and A6 into A1 gives the first result: 


2 2 

( j n 2 -a[v +u -2ucos0 o vcos0cos0-2usin0 o vsin0] 

dvd0d0 = AV COS0 e A7 

This is the equation for the distribution function of a maxwellian 
gas convected with speed u, expressed in spherical coordinates* 

The Explorer 34 plasma probe has a narrow aperature which 
admits only particles for which -A0 < 0 < A0. It is sensitive only 
to the velocity component v cos 0; but [ 0 | < A0 ^ 9° for the 
particles which are detected, so vcos 0 ~ v. The probe integrates 
the velocity distribution over all 0 by virtue of its spin about 
the v z axis , and it integrates over the interval ~A0 < 0< A0 by 
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virtue of its aperature in the 0 direction. Thus, the instrument 
measures 

2 2 

2 -c*(v +v ) 

= Av e x Q ft A8 


where 


Q(0 ,40) = 
o 


40 +2auvsin0 Q sin0 2n + [2auvcos0 o cosO]cos(? 

e cos0 I e d0d0 A9 


<-A0 


The 0 integral in Q(0 ,40) is equal to 


= 2n I o (2o'uvcos0 o cos0) 


A10 


where I Q is the modified Bessel function of order zero. 

Since a is on the order of the reciprocal of the thermal 
speed squared which is generally much less than uv, and since we 
are interested in small 0 q and 0, we can use the asymptotic 
expansion for the B e ssel function to obtain 

2auvcos0 Q cos0 

Q 0 - e 


J 2<juvcos0 o cos0 


AU 
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The remaining integral in (A8) is then 


Q„ = 


V2tt 


0 ^Q'UVCOsPq 


A0 

I 

-A0 


2auv(sin0 o sin0+cos0 Q cos0) 


Vcos© d0 


A12 


We are interested in small 0 (since A0=»9 O ) and small 0 o> (since 
the flow is nearly normal to the spin axis of the detector) , so we can 


expand the trigonometric functions in a Taylor's series to obtain 
A0 


Q 0 = 


- JEL 

J 2 


:auv 


i 


[2o;uv(9 0+1- 

2 


A13 


d0 


2 

Con^leting the square in the exponential and neglecting 0 Q compared 
with 2 gives 


Q 


0 


_ V 2tt 

y2auv 


2a uv 
e 


A0 


-A0 


-2uv(0-0 o ) 2 
e d 0 


A14 


Note that this implies that the 0 dependence of the distribution 

f(v) is approximately Gaussian with a maximum at 0=0 . This is 

o 

a common assumption (Argo et al. 1967). Now we shall distinguish 

two cases. 

1. Interplanetary Medium . In the interplanetary medium essentially 
all of the particles move with (0"0 Q ) small, of order 3 degrees, 
and 0 Q is commonly of the same order, so we can extend the 
limits in (A14) to + ®. 
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In this case 


Q 


0 


2c* uv 

n e 
auv 


A15 


Thus, we find from (A8) and (A14) that in the interplanetary 
medium 


dn _ 
dv 


TT 




2 

-a(v-u) 

e 



‘ 2?T (v-u) 2 
(J) e A16 


where 


■ »/| 

■ / m = .436x10 if cgs units are used. 

2nk 

2. Magnetosheath . If the flow speed is low and the temperature 
is high, as may occur in the magnetosheath, then one cannot 
extend the limits of the integral in (A13). However, the 
integral can be expressed in terms of error functions: 


a) If 0 o < 40 


Q* = 


tt e 


2a uv 


9 2auv 


{erf [,yauv (0 o +o0)] + erf|\/u?uv ( 0 Q -A©) 3 } A17 


and we find that 


dn _ ,436x10 ^ n v 

dV 2 V T ■“ 


- 20 


(er£[.y| (8 o +£,e)] + 


A18 


erfLy /— 2 (0 Q -a0)]3 
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where 


a = S 


u 
k T 


<-u>- 


1.21x10'® X 


v 

u 


A 19 


b) If 9 > A0 

o 


then 


dn _ 
dv 


.-4 


.436x10 N e 
2 


m (v-u) 

2kT i {erf[v^H (0 +a0)]-erfL/£. (0 -A0)]} 


and a is given by (A19). 


A20 
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Figure Captions 

Figure 1. This shows how the telemetry encoding takes place. Data 
for five spins, 1/3 of a spectrum, are transmitted during 
one sequence. The instrument cycles between the proton, 
helium, and fluctuation modes. Three cycles, approximately 
three minutes, are required to obtain a complete spectrum, 
which consists of 15 measurements. 

Figure 2. This shows the definition of a working spectrum and illu- 
strates the fitting procedure. The spectrum shown here (a 
disturbed magnetosheath spectrum) is not typical. It was 
chosen to demonstrate that a general fitting procedure is 
needed because the spectra are not all maxwellian. The 
dashed line between any two bars is obtained by fitting 
three bars to a convected maxwellian distribution; for 
example, the curve between bars 1 and 2 in the working 
spectrum is obtained by fitting bars 1, 2, and 3. 

Figure 3. The upper diagram shows the arrangement of the detector 

relative to the sun sensor in the satellite. The sensitive 
cone of the detector is ~ 2° wide. The rotation direction 
shown here is that which is seen by an observer south of 
the ecliptic plane. The lower diagram shows the sector 
pattern relative to the earth-sun line, which is seen by an 
observer north of the ecliptic plane. The on-board computer 
stores the counts which are recorded as the sensitive cone 
moves through each of these sectors; the counts are recorded 
in the sequence given by the numbering of the sectors. 
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Figure 4. 


Figure 5. 


Figure 6. 


Figure 7. 


Figure 4a shows the ecliptic plane projection of the orbit 
2 trajectory of Explorer 34, and it shows the position of 
apogee for orbits 2 to 20. Figure 4b shows the Y-Z projection 
(projection in a plane perpendicular to the earth-sun line) 
for orbits 2, 6, 10, 14, 18 and 20. 

We have calculated the mean speed <v> and densitity n Q 

which would be observed by the plasma detector if the actual 

distribution function is a convected maxwellian distribution 

— * 

with density n, temperature T, and both velocity u. Here 
we show how <v>/|u| and n Q /n vary as a function of 0 O , where 
0 O is the direction of u as shown in Figure 11. 

A plasmagram, the flow directions, and the magnetic field 
intensity are shown for a four-hour interval during a gep- 
magnetically quiet day (IAGA designation, QQ) . A distinct, 
discontinuous boundary crossing can be seen at 1804 UT., 

There appears to be a very brief excursion into the magneto- 
sheath at 1710 UT. 

The plasmagram above seems to indicate a diffuse transition 
between the magnetosheath and the solar wind; but the flow 
direction and magnetic field intensity indicate that the 
apparent diffuseness is the result of a complex microstructure 
consisting of many discontinuities. The letters S, I, and 
M indicate magnetosheath, interplanetary medium, and 11 mixed' 1 , 
respectively. 
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Figure 10. The solid curves show the theoretical streamlines and shock 
boundary computed by Spreiter et al. (1966) for M» = 6 , 

Y ~ 2 ; the streamlines for y = 5/3 are essentially the same. 
The observed flow direction was into the solid sectors 
which are shown. The sectors shown here have been moved 
slightly in each case so that they fall on the theoretical 
shock boundary. The data were selected so that the observed 
boundary and flow direction were in most cases within 1 Rt* 
of the theoretical boundary. 

Figure 11. This shows the coordinate system which is used for the 
calculations in the Appendix. 
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Figure 8. 


Figure 9. 


An unusual combination of typical transitions between the 
magnetosheath and the solar wind. The discontinuity at 
1409 UT was produced by an interplanetary shock wave that 
presumably drove the bow shock toward the earth past the 
satellite, which was in the magnetosheath just before 
1409 UT. Prior to 1409 UT we see samples of the quiet 
magnetosheath and quiet solar wind; after 1409 UT we see 
samples of the disturbed magnetosheath and the disturbed 
solar wind. The multiple crossings shown here are commonly 
observed; we consider these to be the result of motions of 
a single thin boundary. 

This shows the ratio of the measured mean speed on the 
magnetosheath side of the boundary to the measured mean 
speed on the solar wind side of the boundary, plotted as 
a function of X/D, where X is the coordinate of the satellite 
on the earth-sun line, and D is the theoretical postion of 
the magnetopause boundary at the subsolar point. The data 
were selected so that the interplanetary temperature was 
< 5 x 10^°K; for these data the mean speed is essentially 
equal to the flow speed. The solid curve is the theoretical 
flow speed ratio given by Spreiter et al. (1966). The 
theoretical curve is essentially independent of the Mach 
number and the adiabatic exponent y in the interval shown. 
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